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1 . Introduction 

The problem of stability of the hard-core pinch device has been 
the subject of both theoretical and experimental investigations dur ing 
the recent years (Anderson, Baker et al 1958 ; Anderson, Furth et al 
1 J 58 ; Bickerton and Spalding 1962; Birdsall, Colgate and Furth 1959 , 
i960; Jukes I96I; Reynolds et al 1959 J Tandon and Talwar, 1961, 1962; 
Taylor and Hopgood 1963; and Aitken et al 1964 ). The device consists 
of an annu l a r plasma detached from the central core and the surro unding 
conductirg wall due to opposite axial currents flowing in the core and 
the plasma. The interest in this device stems from the fact that ac- 
cording to ideal (infinite electrical conductivity) hydromagnetie 
theory, the system should be stable against aJ 1 perturbations pro- 
vided the axi al return current through the plasma is less than or 
? 
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equal to the current in the central rod. On the other hand experi- 
ments ( Aitken, 1964) have shown the presence of irreproducible magnetic 
prohe signals indicative of fine scale instabilities in hard core 
geometry where the field distribution satisfying the stability criterion 
based on ideal hydromagnetic magnetic theory were set up. 

There may be several causes for this disagreement between experi- 
ment and theory. The chief reason perhaps is that the idealizations 
de m anded by the hydromagnetic theory are not fulfilled in the experi- 
ment. One, for example, is that the instability may be due to finite 
electrical conductivity of the plasma which allows a relative slip- 
page between the lines of force and the fluid, thus removing the con- 
straint imposed by infinite conductivity. This possibility has 
recently been explored by various investigators (Bickerton and Spalding, 
Tandon and Talwar, Taylor and Hopgood) who have shown that in the limit 
of vanishingly small electrical conducivity, the hard core pinch 
model is, in general, unstable. This gives one the feeling that the 
observed instability is due to finite conductivity, although a final 
answer can only be given when calculations are done with finite (and 
not vanishingly small) conductivity. 

The purpose of the present paper is to investigate the instability 
of an idealized hard core pinch configuration, regarding the plasma 
pressure to be anisotropic rather than scalar. This happens in the 
case of a dilute plasma subject to a high magnetic field. In the 
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simplest situation the plasma pressure is different in at least two 
directions viz parallel and perpendicular to the magnetic field. It 
is well-known that in case the plasma pressure is initially (before 
perturbation) isotropic in a static plasma, the MED approximation gives 
the most pessimistic result regarding stability, i.e., if we obtain 
stability in the MHD approximation, the stability is certa in ly guaran- 
teed in the more exact kinetic theory in the (M/e) -» o limit as well 
as in the Chew, Goldberger and Low approximation (1956). If, on the 
other hand, the equi3J.br ium state is characterized by an anisotropic 
pressure, new types of instability may arise - a typical example being 
the "mirror" and the "fire-hose" instability for plane waves in a 
homogenous plasma carrying a uniform magnetic field. 


2. Formulation of the Problem 



Consider a low density plasma contained between the two 

coaxial cylinders formed hy the metallic core of radius A 0 R 

and a regid perfectly conducting wall of radius A 2 R. Axial 

currents of strength J and J are respectively flowing in opposite 

c p 

direction through the central conductor and the plasma as a result of 
which the plasma will detach itself from the conducting central rod 
and the wall forming an equilibrium model of hard-core pinch as shown 
in Fig. 1. Let R and A x R respective denote the radii of the inner 
and outer plasma-vacuum interfaces A 0 , and -.A 2 axe constants 

such that 

A 0 < J } 


we further assume that there is a strong axial, magnetic field B^ 
prevalent inside the plasma mating the plasma pressure anisotropic. 
Let B 2 and B 3 be the axial magnetic fields in the inner and the outer 
vacuum regions respectively. 

Assuming that there are only surface currents on the two plasma- 
vacuum interface, we can write for the equilibrium fields as 



A 0 £0 

< A 1* 


(2) 
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We may egress the magnetic fields in various regions in terms 
of the toroidal field B. at the inner (r = E) plasma -vacuum interface. 
We thus have. 




( 3 ) 


and 


a 



J 

- A . 
J 1 




here B & denotes the toriodal field at the outer plasma-vacuum inter- 
face. In equilibrium the total stress (plasma plus magnetic pressure) 
should he continuous at each hounding surface. This gives 



ok 


(5) 


and 
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oir ^ ^ R. 


(6) 



It can readily "be seen from equations (5) and. ( 6 ) that the 
field parameters in equilibrium should satisfy the relation 


O. V?" ( 7 ) 

1 +^ cat \ 

Further, the equations governing the motion of the non-dissipative 
p lasma are written as 

As . . „ ( 8 ) 

$ - ~ v i* — 

_ q. (f v) (9) 

!i . v* W* &) (10) 

H t “ 

^ - 41 ? & (ll) 

r-T n (12) 

V. 2 - Q 


(p^ vv, B and respectively denote the density, velocity magnetic 
field and current density vectors at a point), p the pressure tensor. 
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which is diagonal in a local rectangular coord i nate system with 
one axis along the magnetic field and invariant under rotation about 
B is given by 

= \>„ th-* -v h.Q, - ' 1 * («) 

Here n is a unit vector pointing along the magnetic field and 1 
signifies the unit dyadic p xx and p^ represent the gas pressures 
along and perpendicular to the fiels respectively. 

Neglecting the heat flow tensor, we further have the double 
adiabatic equations (Chew, Goldberger and Low 1956) 



( 14 ) 


( 15 ) 


It can readily be scan from equation (8) that in the initial 
state p xx and p^ are constants . 

In the two vacuum regions, the equations to be satisfied are 
(ll) and (12) with _ 0 

We study the stability of the static configuration characterized 
by the above equilibrium parameters and as depicted in Fig. 1. by 
imparting an infinitely small velocity perturbation v to the system. 
The various resulting perturbations, like the change in gas pressure. 
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magnetic field, current density etc. are assumed to "be of the form, 

(function of r) x e jcp. { nt + i (m 9 + k z) ] (l6) 

where n is the growth rate parameter, m the azimuthal number and k 
the wave number of the perturbation in the z-direction. A dispersion 
relation relating n, the parameter determining the stability, with 
other physical qualities is then obtained by adopting the us ual nor mal 

i 

mode analysis .with suitable boundary conditions. 
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3- Equations Governing the Perturbed State and Their Solutions 

With the perturbations of the form (l6) the linearized pertur- 
bation equations for the plasma can be written as 

s > (17) 

(18) 

'W.Sfc - C h - v ) - - hV')L (19) 

■= <0 ( 20 ) 

V.'bb - o (21) 

^ - V*y.^ (22) 

* % - K - %v --) (25) 
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where 


2 * , s f, s 3 


and U denote 


the perturbations^ of first order of smallness^ in various physical 
parameters. 

The perturbed equations for the vacuum regions are 


x?y.SB ce \ 


= ° 


where b 6 tf> represents the change in the magnetic field in the two 


vacuum regions. 


(a) Solution for plasma region 
Writing 


$* 1 . s>= 


we obtain from equations (17) to ( 23 ) that 


U = b, ,^4 e 
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where 


( 29 ) 


A - 




(30) 


and 4 . satisfies the relation 
J Z 




( 31 ) 



Solution of equation (30) is of the form 
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V- + 


where A x and Aa are arbitrary constants. 

(b) Solution for vacuum regions 

The change in magnetic field, ^ inside the vacuum 

regions is given by equations (24) and (25). Thus sj' #) 

is 

derivable from a scalar potential X. , satisfying the relation 


2 § ^ , v% 


This equation when combined with equation ( 25 ) leads to 




/5o 


X = [c,I^ O'*-) + + ' (r>e 


Since there are two vacuum regions, separate boundary condi- 
tions are required for the two regions. The boundary conditions for 
the inner vacuum region are 


N. § - o 


Ar S A, 


= 


(d: a,- ^4- i 
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where M is the normal vector to the surface and -5 is the dis- 
placement at the inner-plasma vacuum interface. 

Similarly the boundary conditions for the outer vacuum regions 
can be written as 


N . 5 — o 


oJc~ nj - 


(58) 


and 


Im] - 


oJtr 


A/ 


- A . R 


3T 


it? 


where 



represents the displacement at the outer plasma vacuum 


interface . 

With the help of boundary conditions (57) ani (58), we write 
expressions for the change in the magnetic fields at the two plasma- 
vacuum interfaces as 




(*»*•*■)£ 
S x ) Qn> (.*< I 'l’0 J 


( 59 ) 


and 



where 


Qnv {h%) 


C (^) iL (?) 


(4l) 


a. 


X SI fe.& 


( 42 ) 


Suffix I and II represent the parameters for the inner and outer 
plasma-vacuum interfaces respectively. 
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4. Dispersion Relation 

In order to obtain the dispersion relation the pressure ( kin etic 
+ magnetic) balance condition must be satisfied at the two plasma - 
vacuum interfaces in the perturbed state. These conditions in the 
linearized form can be written as 


'h 






e> x . 5E 1 


4ir 




eJt ^ 


and 


l h 


b,.L 6 







TT 


ok n.^ 


Conditions (43) with the help of equations (2) - (4), (27) - (33) 
and. (39) to (4l) give the following characteristic equation relating 
n 2 with the various equilibrium parameters 


F +- F 4- - o 


(44) 


where 
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£ ^ ^ ^ ^ C/M) ~ C*>) ^ ( A W ) 

CW K~ CM) - C ft) £~ CM) 


(50) 


P C.1 m) - 3Tm c^) Ky " Cm) - ^ c^) c a ^) 

x^c^) ^ Cm) — ^*»(m) 


(51) 


‘jj.s j^ft. a^l "the functions Q 


■with 


are given by the relation (4l) . 
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5 . Neutral Stability Curves 

The dispersion equation (44) has to he solved, for the parameter 
n 2 in order to decide whether the configuration is stable or not. 

For any set of given physical parameters we have to incorporate the 
equation (45) with equation (44) while eval ua ting n 2 . The equations (44) 
and (45) are extremely unwieldy and we shall he content with investigat- 
ing the question of marginal stability (n = o) for anisymmetric per- 
turbations (m = o). The eigenvalues n 2 for a non-dissipative, static 
plasma configuration (as is under investigation) are known to be real 
(positive or negative) so that the neutral stability curves define 
the transition from stable to unstable regime. 

Numerical calculations for the roots F of the quadratic equa- 
tion (44) were done on computer taking A = 0.5, A, = 2, A = 4 
and b 2 = b x , b 3 = Jl - a^. With a 2 as positive (equation 52 ) 
one root of the quadratic equation was plotted against x (=kR) to 
obtain neutral stability curves in a variety of cases, (figures ( 2 ) - 
( 5 )). each figure curves (a), (b), (c) are for a = 0.2, 0.4 and 
0.6 respectively. Figures ( 2 ) and ( 3 ) are for b 2 = b x = o and a 2 = 0.1 
and 1.0 whereas b 2 = b x = 1, a 2 = 0.1, and 1.0 are the respective 
values for figures (4) and . ( 5 ). Results were also obtained regarding 
a 2 in equation (32) as negative and are plotted in figures (6) - ( 9 ) 
t aking b 2 ^ b x = o, and b 2 = b x = 1 for same set of values for | a | 


and a 2 . 
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ABSTRACT 


The stability of a hard-core pinch model, characterized by an 
anisotropic pressure, is investigated using the usual normal mode 
technique. The dispersion formula is obtained with the "double- 
adiabatic" hydr omagne t i c equations and discussed for marginal 
stability for axisymmetric perturbations. 
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CAPTIONS TO FIGURES 


Fig. 1 
Fig. 2 

Fig. 5 

Fig. 4 

Fig. 5 

Fig. 6 

Fig. 7 


Equilibrium model of hard-core pinch. 

Neutral stability curves for b^ = b^ «= 0, and b^ 2 = 0*9 , 
and A = 0.5, A, = 2, A =4. Curve a, b, c respectively 
correspond to a = 0.2, 0.4, and 0.6. 

Neutral stability curves for b =b n =b =0, A = 0.5, 

A^ = 2, and A^ = 4. Curves a, b, c respectively correspond 
to a = 0.2, 0.4, and 0.6. 

Neutral stability curves for b =b, =1, b 2 =1.9 and 

2 1 3 

A = 0.5, A. = 2, A =4. Curves a, b, c respectively 

O X 2 

correspond to a = 0.2, 0.4, and 0.6. 

Neutral stability curves for b^ = b^ = b =1, and A q = 0-5, 
A-, = 2, A =4. Curves a, b, c respectively correspond to 
a = 0.2, 0.4, and 0.6. 

Neutral stability curves for b = b-= 0, b 2 = 0-9 with 

2 l 3 

A q = 0.5, A^ = 2, A^ = 4. Curves a, b, c respectively 
correspond to a = .2i, .4i, 0.6i. 

Neutral stability curves for b = b n = b =0 with A = 0.5, 

2 1 3 o 

A n = 2, A =4. Curves a, b, c respectively correspond to 

x 2 

a = .2i, 0.4i, 0.6i. 

Neutral stability curves for b = b n = 1, b 2 = 1.9 with 

2 J- 3 

A q = 0.5, A-^ = 2, A^ = 4. Curves a, b, c respectively 
correspond to a = .2i, 0.4i, 0.6i. 


Fig. 8 
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Fig. 9 Neutral stability curves for b = b. = b =1, with A =0.5. 

2 1 3 o * 

= 2 , A^ = 4. Curves a, b, c respectively correspond to 
a = 0.2i, 0.4i, 0.6i. 
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Figure 3 





















